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Torrefied biomass provides opportunities for an alternative, renewable fuel in the energy market. Much 
work on torrefaction is based on weight loss transients without much insight into the evolution of volatile 
species. A previous model from our group has been developed that predicts the degradation of raw 
biomass into specific chemical species during torrefaction (water, acetic acid, carbon monoxide, carbon 
dioxide, formic acid and furfural) in three consecutive reaction steps. The goal of this work is to understand 
how the kinetics of torrefaction change with temperature. Kinetic parameters for aspen wood were 
obtained for the first 90 min of 260-300 °C (1000 °C/s) torrefaction. The pre-exponential factors were 
3.32E9, 1.43E11, and 2.08E14mirr\ and the activation energies 1.05E2, 1.27E2, and 1.72E2kJ/mol for 
the three reactions. Kinetic parameters found correspond well to similar values from global weight loss 
studies, and suggest that increased torrefaction severity leads to progressively more recalcitrant forms 
of torrefied biomass. These continuous production traces for the volatile organic species yield insight 
into torrefaction, and require careful consideration of torrefaction time and temperature to design the 
desired solid fuel product. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

Currently 42% of the electricity generated in the United States 
comes from coal combustion [1,2]. Coal is known to produce one- 
third of the C0 2 emissions in the United States and releases other 
harmful pollutants such as CO, SO*, NO x , particulates, and mercury. 
Because of the growing concern of global warming, the need to find 
renewable C0 2 neutral feedstock is immediate. Biocoal (or torrefied 
biomass) [3-6] can be an immediate solution in coal utility boilers. 
This renewable fuel can be used as a drop-in replacement for coal 
firing due to its high energy density, similar properties to coal, low 
emissions and potential for C0 2 neutrality [3,6,7]. 

Biocoal is produced from the torrefaction of biomass, at temper¬ 
atures from 200 to 300 °C in an oxygen free environment. During 
torrefaction, mostly the hemicellulose portion of the biomass 
degrades through a complex coupled chemical-kinetic-heat-and- 
mass-transfer process 7-9]. There are small contributions from 
cellulose and lignin during this process; with the contribution 
of cellulose increasing significantly at temperatures above 300 °C 
[8,10]. When the biomass is heated up quickly in an inert atmo¬ 
sphere, biomass is dehydrated, the viscoelastic properties are 
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changed, and the calorific value increases 6,11]. In addition, the 
energy requirements to reduce the size biocoal is substantially 
smaller, up to an order of magnitude, when compared to raw 
biomass [4,5,7,12,13]. Biocoal is then able to be a drop-in replace¬ 
ment for pulverized coal. 

In order to have a better control of the torrefaction process 
and develop more efficient reactors for biocoal production, the 
torrefaction of biomass and the kinetics of species evolution, and 
changes in the residual biomass need to be better understood. 
Much gravimetric work has been done to model changes during 
torrefaction such as parallel reaction [7,14-16] and consecutive 
reactions models [10,17-20]. These models are useful for targeting 
throughput in equipment design, but provide little insights into the 
value-added benefits of torrefaction such as the improved fuel and 
physical properties over raw biomass mentioned above. A detailed 
understanding of the feedstock degradation during the torrefaction 
process, or the gas species evolution, is necessary for prediction of 
product quality changes with reaction conditions. 

One specific approach to modeling torrefaction kinetics is to 
develop equations around weight loss data. One such model that 
is well accepted in literature is the two stage degradation of hemi¬ 
cellulose (xylans) developed by Di Blasi and Lanzetta [17]. This two 
stage model first degrades to a solid intermediate through rate con¬ 
stant /<i, and volatiles through kv \, followed by a sequential step to 
torrefied biomass through k 2 and more volatile products through 
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kv 2 . This model predicts one stage of degradation at early times, 
and one stage for degradation over long processing times. The rate 
constants for each stage are also combined to give a global rate 
constant for that step ( k\ = fc B + kv \, k 2 = k c + kv 2 ). This was further 
developed by Di Blasi and Branca in subsequent years to include 
a similar third degradation stage (constants now /< B and kv lt /< D 
and kv 2 , and kc and kv 3 , respectively) when whole wood species 
(beech) were studied instead of xylans [21,22]. Further, these rate 
constants were combined to predict the pseudo-reactions where 
/<i is the sum of the constants for volatile production and the pro¬ 
duction of the solid intermediate ( k\ = kv\ + /< B ) in reaction one. The 
six rate coefficients were therefore combined and analyzed as three 
global constants (/<i = kv i + /< B , k 2 = kv 2 + /< D , k 3 = kv 3 + k c ) from 528 K 
to 584 K (6000 s), and 573 K to 7081< (up to 2000 s). The first reaction 
(/<i) was reported to dominate from 528 K to 593 K, the second (k 2 ) 
from 528 K to 708 K, and the third (k 3 ) from 603 K to 708 K. While 
these models may be convenient for estimating material through¬ 
put in equipment design, they do not provide any process insight 
into product quality, or any information on a chemical level. 

Several studies have investigated modeling of torrefac- 
tion and coupling information about volatile species evolution 
[ 10,18-20,23,24]. Prins et al. studied the torrefaction of willow over 
225-300 °C using the Di Blasi and Lanzetta model discussed above 
[10]. This was quickly followed by a subsequent study investigat¬ 
ing the product distribution that is observed [23]. These were later 
combined to accurately predict the volatile product evolution over 
a range of temperatures and times [19]. More recently, Nocquet 
et al., took a similar approach using beechwood, cellulose, lignin, 
and xylan standards [20,24]. In the kinetic study, the two-stage 
Di Blasi and Lanzetta model was used to model the decompo¬ 
sition of the cellulose, hemicellulose, and lignin fractions of the 
biomass, with a stoichiometry parameter for the production of the 
respective volitile species [20]. With a correction facor to account 
for the interaction of the related consituents of virgin wood, they 
were able to accurately predict the volitile evolution at descite 
points, and also match a continuious mass loss signal. In order to 
move beyond a descrete point understanding of volitile torrefac¬ 
tion species, however, our group previously investigated a model 
based strictly around gas-phase species kinetics [18]. 

The previous study predicts the first 90 min of torrefaction at 
300 °C, and the yield of products in a three-step torrefaction model 
[18]. This model was able to accurately predict the amount of var¬ 
ious chemical species evolved during torrefaction and was based 
on the degradation of the hemicellulose material at 300 °C. The 
reaction model is shown below. 


+ Y) a i p i 

i 

(i) 

rYr 2 + Y^ p ‘ 

i 

(2) 

R 2 \c + Yn p i 

(3) 


In this model the initial hemicellulose (H) is degraded with reac¬ 
tion rate constant k\ into a reaction intermediate R\ and forms a 
suite of products {P/,i}. The first reaction intermediate R i is then 
degraded with rate constant k 2 into a second reaction intermediate 
R 2 and produces more of these products {P; )2 }. The second reaction 
intermediate is then degraded with rate constant k 3 into the biocoal 
product (C) along with an additional release of reaction products 
{P l 3 }. This model accounts for the production of major torrefaction 
products through characteristic mass spectrometer (MS) ion frag¬ 
ments of: water (18 m/z), carbon monoxide (28 m/z), formic acid 


(29 m/z), carbon dioxide (44 m/z), acetic acid (60 m/z) and furfural 
(95 m/z). 

This paper aims to deepen this model by applying this model 
over a range of torrefaction temperatures, and investigate the kinet¬ 
ics of torrefaction purely on the gas-phase reactions. In this work 
we take a lumped parameter assumption, where that the reaction 
rate coefficients are the same for each compound evolved, but the 
intensities of the specific species change with a stoichiometric coef¬ 
ficient {a,p,y). These coefficients represent a relative amount of 
volatile product generated from each consecutive reaction stage, 
and are unique for each species—each of the six compounds inves¬ 
tigated have a stoicheometeric parameter for each reaction. If this 
model is valid over the temperature range, then the stoichiomet¬ 
ric coefficients should be independent of temperature for each of 
Eqs. (l)-(3). This assumption was validated in this research. This 
approach to predicting torrefaction is novel, but was previously 
investigated at only at 300°C [18]. By obtaining gas-phase kinet¬ 
ics, accurate predictions of the volatile yields can be obtained for 
both isothermal and non-isothermal conditions. The dynamic evo¬ 
lution of these species was investigated, and unique implications 
of interpreting torrefaction through weight loss data in the absence 
of these novel results is discussed. 

2. Materials and methods 

As in our previous study 18] we used aspen biomass for this 
current research. Aspen samples were debarked then dried in an 
oven at 105 °C until no difference in mass was observed. A Thomas 
Wiley® knife mill (NR. 3557524 359264) was used to reduce the 
size of the wood particles. To obtain particle sizes of 0.5-0.6 mm a 
W.S. Tyler Rotap (model RX-29, serial 9774) was used. The aspen 
samples were roughly cylindrical with diameters of approximately 
0.5 mm and length ranging from 1 to 3 mm. 

2.1. Torrefaction experiments 

To obtain kinetic data, the dried aspen wood samples were 
torrefied using a CDS Analytical micro-pyrolysis unit (model 
5200HP Pyroprobe). The aspen was torrefied at temperatures of 
260 °C to 300 °C at 10°C intervals for 90 min. The volatile species 
evolved were continuously analyzed using a Trace GC Ultra (Model 
K8880181) gas chromatograph by ThermoFisher and a quadrupole 
mass spectrometer (Trace DSQII, ThermoFisher). The GC was fit 
with an inert, non-retentive, 10 m fused silica column (Restek, Rxi 
Guard Column). This eliminated any interaction with the com¬ 
pounds in the GC, and acted as a transfer line from the torrefaction 
reactor to the MS. 

Aspen samples of 0.1-0.2 mg were weighed to within an accu¬ 
racy of 1 jjug using a microbalance (Citizen Scales Inc., Model 
CM5). The samples were then placed in a quartz reaction tube 
(25 mm x 2 mm OD) and held in place by two small pieces of glass 
wool on either end of the reaction tube. The samples were loaded 
into the micro-pyrolysis unit such that the wood chips were within 
the center of the heating zone and rapidly heated to the desired 
temperature. Rapid heating (1000°C/s) was accomplished using 
a platinum heating filament that surrounded the quartz reaction 
tube. Ultra high purity helium (99.999%) was used as the carrier 
gas in order to create an inert atmosphere within the reaction zone 
and to transport the gaseous and vapor products into the inlet port 
of the GC through a heated stainless steel transfer line (300 °C). The 
GC column was kept at 300 °C for all experiments and the flow of 
helium gas kept at 1.5ml/min (superficial velocity of 50.9 cm/s). 
After the samples passed through the GC, they then entered the 
quadrupole mass spectrometer where the ion fragments’ intensi¬ 
ties with respect to time were recorded. The mass spectrometer 
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measured the continuous signal of ions from 15 to 350 m\z for 
90 min. Data were then exported to a Microsoft Excel sheet for 
analysis. The measured data were normalized by dividing the MS 
intensity measured at each time point by the initial mass of aspen. 

2 . 2 . Model fitting 

As mentioned, the proposed model (Eqs. (1 )-(3 )) is a consecutive 
3-step degradation model for the early stages of torrefaction (up 
to 90 min). This model was developed for the main products from 
torrefaction: CO, C0 2 , H 2 0, organic acids (acetic acid and formic 
acid), and furfural [9,18,23] (although some studies report that 
methanol, methane, formaldehyde, and hydrogen are also observed 
[9,23-29], they were excluded here because they were not formed 
in any detectable level in our studies). Measuring these compounds 
over the course of torrefaction can yield information about the gas 
phase formation kinetics. In addition, changes in recalcitrance of 
the torrefied biomass, and changes in elemental composition can 
be inferred. 

From our previous study [18] it was shown that at early times 
the production rate from Eq. (1) is dominant, while at later times 
Eq. (3) is the most important for product formation. An estimate for 
the parameters in Eq. (1) ( H 0 (Xi , k\) during early time (0.5-2.0 min) 
and the parameters for Eq. (3) (H 0 y it k 3 ) at late times (>60min) 
were obtained by adjusting the parameters to match the pro¬ 
duction intensity (stoichiometric coefficients) and dynamics (rate 
constants) over the specified times. The final parameters for Eq. (2) 
(H 0 /3j, k 2 ) were then fit using a nonlinear, generalized reduced gra¬ 
dient method. The fitting routine was performed in Microsoft Excel 
around reducing the sum of squared error between the model and 
all measured data. 

This method assumes that torrefaction is a set of devolatilization 
reactions that take place over a range of temperatures. Further, this 
implies that the stoichiometric parameters (a, y) are indepen¬ 
dent of temperature. In addition, the model predicts lumped groups 
of product formation (i.e. all species have the same rate constant 
for each reaction step). These assumptions in the formation of the 
model are important and yield constraints about the system, and 
chemical significance. Because the same reactions are assumed to 
occur over the temperature range, a constraint is imposed that the 
stoichiometric parameters should not deviate with temperature. 
Validation of this is investigated in a following section. Further, 
the kinetics of torrefaction have been shown to closely follow tra¬ 
ditional Arrhenius behavior [10,11,15-17,22,28,30,31] and impart 
further insight into kinetic aspects of the modeling. Initial estimates 
for the rate constants described above were used to obtain an initial 
set of Arrhenius parameters. These values were then used to obtain 
a new set of rate constants that are strictly linear on a In k vs. 1 IT 
basis. These values were then allowed to change along with the stoi¬ 
chiometric parameters (a, y) in the iterative process of reducing 
the sum of squared error between the model and the experimental 
data. This method was followed for all species and temperatures. 

All six volatile species were examined over the temperature 
range of 260-300 °C, however, when the torrefaction temperature 
decreased to 260 °C, carbon monoxide was not generated in suffi¬ 
cient quantity to produce a clear signal in the detector. A very low 
signal to noise ratio (S/N) makes it impossible to distinguish from 
the background noise, and these data were therefore excluded from 
the set. 

3. Results and discussion 

3.1. Experimental and modeling results 

Fig. 1 shows typical traces of all compounds at 20 °C inter¬ 
vals (to simplify the presentation) to demonstrate the different 
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Fig. 1. Experimental transient trends for 1 8 m\z (H 2 0 ), 28 m\z (CO), 29 m\z (HCOOH), 
44 m\z (CO2), and 60 m\z (CH3COOH), and 95 m\z (OC4H3CHO) at varying tempera¬ 
tures. 


behaviors observed in the formation of gaseous and vapor products. 
Remaining data traces are shown in the supplemental materials. 
These traces are directly related to the rates at which these compo¬ 
nents are produced. The figure also shows that as the temperature 
increases, the characteristic reaction times decrease. At early times, 
during the initial decay, we can observe the highest intensities 
for the products at the highest temperatures. As the temperatures 
decrease these production rates slow, and the production is pro¬ 
longed over later times. This behavior flattens out the production 
curves, and thus demands much longer times to reach the same 
production amounts. From the figure, the steady decrease in pro¬ 
duction rate with temperature is intuitive and suggests simple 
temperature dependence. The production rate data decrease sig¬ 
nificantly over the 40 °C investigation range. This suggests that 
the process is very sensitive to both temperature and time, and 
any relative decrease in temperature would require a much longer 
processing time to reach the same feedstock degradation. 

Except for very early times, all compounds appear as a mono¬ 
tonic decay that decreases with decreasing temperature, but 
require several reactions to describe both the early and late time 
production rates. The acetic acid and furfural behavior is more 
complicated, and exhibits a slightly delayed production. These 
compounds are dominated by a delayed apex production peak 
shortly after the initial rate. The modeling of these traces, as shown 
in Fig. 2, demonstrates that the initial monotonic production (Eq. 
(1)) of acetic acid and furfural is quickly overtaken by a larger 
production in the subsequent reactions (Eqs. (2) and (3)). This 
generates a production apex after the initial time. Like the other 
traces, this peak decreases in size with decreasing temperature as 
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Fig. 2. Model fitting for all compounds at different torrefaction temperatures. The 
experimental data are shown as solid lines. The model at each temperature is shown 
as dashed lines. 


the production delayed in time enough that it is masked by the 
slow kinetics. Fig. 2 shows typical model fits for all six compounds 
at 20 °C intervals. Good agreement between the model and the 
recorded data was observed over the temperature range investi¬ 
gated. Remaining transient sets are displayed in the supplemental 
section. 

3.2. Model validation 

It is important to note, again, that the model was assumed to 
apply over the whole temperature range—meaning the stoichio¬ 
metric coefficients (a, p, y) should be independent of temperature. 
In addition, the rate coefficients are species independent - mean¬ 
ing the model predicts lumped groups of product formation. The 
kinetic parameters were fit to the experimental data as explained 
above and for a , p, and y- To examine the variability in these 
stoichiometry parameters, the average value of ot determined for 
all species and temperatures, and a/a aV g was plotted as a func¬ 
tion of temperature. These results are shown in Fig. 3. The figure 
shows the scatter in values for a away from the average. Although 
there is some variability, the value is approximately constant with 
temperature. Similarly the contribution from (Eq. (2)) and (Eq. 
(3)) were investigated as functions of temperature, and in rela¬ 
tion to the production from (Eq. (1 )) through a. Fig. 4 shows values 
for plot and y\ot over the experimental temperature range, again 
supporting that the values are temperature independent. These 
confirming aspects of this work can lead to further chemical insight 
and development from a lumped product model to a more detailed 
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Fig. 3. Values of alpha for each component over 260-300 °C. 


P,260-300C 



Torrefaction Temperature (°C) 


y, 260-300C 



X 60 

• 95 

■ 28 

A 29 

X 44 

♦ 18 


Torrefaction Temperature (°C) 
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reaction degradation mechanism. In conclusion, all values obtained 
for a/aavg, PI & and y\ot validated that the stoichiometric coeffi¬ 
cients are temperature independent, because we are dealing with 
the same reactions. 
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Table 1 

Kinetic parameters for various feedstocks in a two-stage degradation scheme. 


Feedstock 

Al (min -1 ) 

i A2(min -1 ) 

i El (kj/mol) 

E2 (kj/mol) 

Source 

Xylan 

2.17E7 

2.30E4 

76.6 

54.8 

[17] 

Beech 

1.6 E6 

3.6 Ell 

76 

143 

[27] 

Straw 

1.46E6 

3.26E3 

64.6 

47.3 

[30] 

Corn stalks 

3.78E8 

1.62E5 

91.5 

65.2 

[30] 

Willow 

1.49E6 

6.6E11 

76.0 

151.7 

[10,19] 

Wheat straw 

2.09E6 

2.60E5 

71.0 

76.6 

[24] 


3.3. Arrhenius parameters 

Fig. 5 shows the relationship between the three reaction rate 
constants (as obtained from our fitting procedure) and temper¬ 
ature. The error bars represent the standard deviation between 
the specific rate constant for all species at the specified temper¬ 
ature, showing very good agreement (with R 2 for the three trends 
0.99, 0.99, and 0.98, respectively). These data collected from the 
chemical species model employed in this work show, as others 
have with only weight loss curves, that the kinetics of torrefaction 
closely follow Arrhenius behavior [10,11,15-17,22,28,30,31]. The 
pre-exponential factors are 3.32E9, 1.43E11, and 2.08E14min -1 , 
respectively. The activation energies are 1.05E2, 1.27E2, and 
1.72E2 kj/mol, respectively. This Arrhenius behavior demonstrates 
that the three reaction model is able to capture the temperature 
dynamics quite well. These kinetic parameters can be used to pre¬ 
dict changes in rates at reaction conditions and hence control the 
torrefaction process and its product distribution. 

The proposed model assumes that each reaction produces a 
group of products through a lumped kinetic rate constant, or, 
the unique reaction rate constant is responsible for production of 
all species in that reaction step. The determined kinetic param¬ 
eters, then, are related to an equivalent scheme and parameters 
based on weight loss measurements. The two-stage Di Blasi and 
Lanzetta model [17] has been applied to different types of feedstock 
include xylans [17,20,30], woody feedstock [10,19,20,22,30,32], 
and straws [28,33]. Table 1 shows a summary of kinetic param¬ 
eters for the various feedstock. The activation energies for the 
three stage Branca and Di Blasi model [21,22] were 76.2, 142.8, 
and 43.8 kj mol -1 , respectively, and the pre-exponential factors 
were 2.69E4, 5.91 E9, and 10.0 s- 1 (1.6E6,3.6E11, and 6.0E2min- 1 ) 
respectively. While the present study predicts higher activation 
energies are required for decomposition and a progressively more 
recalcitrant solid structure as torrefaction progresses, the overall 
magnitude of the parameters are similar, and should differ with 
biomass source. 

Over the temperatures and duration of torrefaction studied, 
the present work has reasonably predicted the product forma¬ 
tion of the six major gas/vapor products. This accounts for any 
changes, by difference, in the torrefied biomass. In addition to 


Table 2 

Relative stoichiometric coefficients produced from reaction one (a), two (p), and 
three (y) ±one standard deviation. 


Compound 

a, Unity 

P, Norm 

y, Norm 

H 2 0 

1.0 ±0.2 

1.6 ±0.1 

0.1 ±0.0 

CO 

1.0 ±0.3 

3.6 ±0.1 

4.4 ±0.1 

co 2 

1.0±0.1 

1.7 ± 0.1 

4.5 ±0.1 

HCOOH 

t.0±0.1 

3.4 ±0.2 

4.7 ±0.3 

CH3COOH 

1.0±0.1 

12.0±0.5 

14.0 ±0.7 

OC4H3CHO 

1.0 ±0.2 

5.9 ±0.3 

6.0 ±0.4 


accurately predict the observed experimental data, the predicted 
kinetics are repeatable and closely follow traditional Arrhenius 
behavior. The observed energies have a similar magnitude to val¬ 
ues that others have obtained with similar weight-loss schemes 
discussed above. The activation energies proposed in this work 
predict that at early times the biomass structure is relatively eas¬ 
ily broken down compared to Eqs. (2) and (3). As torrefaction 
progresses, however, the subsequent reaction steps require more 
energy. This increase from Eq. (1) to Eq. (2), and from Eq. (2) to 
Eq. (3) implies that the solid intermediates are becoming pro¬ 
gressively more recalcitrant with torrefaction severity. Eleat and 
mass transfer effects have been shown to be negligible in previ¬ 
ous work [ 18 ]. Therefore, the instantaneous heating rates employed 
are able to decrease any impacts or complications of heat ramping 
involved in prior TGA methods. Finally, the modeling of experimen¬ 
tal data supports our assertion that the same global reactions are 
taking place over the range of torrefaction temperatures—or that 
the stoichiometric ratios of product formation are independent of 
temperature. 

3.4. Reaction contributions in product profiles 

The stoichiometric coefficients obtained from our work provide 
a way to compare the production contributions from each reac¬ 
tion step. Note that the coefficients in the above mechanism are 
specific to each component, but independent of temperature. To 
more readily compare these coefficients, they are normalized to 
the amount of product formed from Eq. (1). These relative con¬ 
tributions are shown in Table 2. As can be seen from the table, 
the majority of water is produced from the first (a = l) and sec¬ 
ond {fila = 1.6) reactions. Very little production {yla = 0A) from 
the third reaction takes place at long times. In stark contrast to 
that, most compounds are only produced in significant quanti¬ 
ties at later reaction times; for example acetic acid has small 
contribution from Eq. (1) (a = l), with a major contribution from 
Eq. (2) {fi/ot = 12.0) and Eq. (3) {yla = 14.0). This suggests that 
very careful consideration of torrefaction severity and duration 
is required to match the desired process outcomes with process 
conditions. 

If these results are examined through the lens of an anticipated 
product distribution, the relationship between weight loss, product 
evolution, and product characteristics can be realized. It is known 
that the gas phase products are comprised mostly of water, acetic 
acid, and C0 2 (for example Ref. [9] indicated for willow the prod¬ 
ucts are 89.3% water, 4.8% acetic acid, and 3.3% C0 2 ; other literature 
values are in the same magnitude but differ slightly with biomass 
source). Using our stoichiometric parameters and the product dis¬ 
tribution of Ref. [9], the product evolution profiles were found 
and are shown in Fig. 6. If total weight loss is examined alone, it 
would appear that only a few minutes are required to complete the 
majority of the torrefaction process. However, an examination of 
the chemical species evolution reveals that long processing times 
are required to evolve carboxylic and other group compounds. 
Removal of these groups from the raw biomass is required for 
the biocoal to resemble the characteristics of fossil coal [3-7,34]. 
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Fig. 6. Total product profiles for water, acetic acid and carbon dioxide with a product 
distribution. 

Similar delayed feedstock decomposition and gaseous/vapor evo¬ 
lution profiles have been previously observed, but not modeled or 
extensively characterized [ 10,19,22,35]. It is our assertion here, and 
from our previous study [18], that while these late-time degrada¬ 
tion steps may appear negligible from a weight loss viewpoint, they 
are an important consideration consisting of removing primary 
organics. These longer reaction times lead to preferable changes 
in torrefied biomass properties. 

4. Conclusions 

This study demonstrates that the three step model can accu¬ 
rately predict gas-phase product kinetics over the initial stages 
of torrefaction (up to 30% mass loss) from 260 to 300 °C. 
The pre-exponential factors obtained were 3.32E9, 1.43E11, and 
2.08E14min _1 , and the activation energies were 1.05E2, 1.27E2, 
and 1.72E2kJ/mol for the three reactions, respectively. At early 
times, the evolution of chemically formed water dominates, while 
later times are required to remove the organic compounds such 
as acetic acid. Understanding torrefaction product evolution may 
lead to better prediction of product quality from varying or non- 
isothermal processing conditions, a topic to be investigated in 
future work. 
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